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Abstract 



The initial- and boundary- value problem for the Benjamin-Bona-Mahony (BBM) 
equation is studied in this paper. The goal is to understand the periodic behavior 
(termed as eventual periodicity) of its solutions corresponding to periodic bound- 
ary condition or periodic forcing. To this aim, we derive a new formula represent- 
ing solutions of this initial- and boundary- value problem by inverting the operator 
d t + ad x — jd xxt defined in the space-time quarter plane. The eventual periodicity of 
the linearized BBM equation with periodic boundary data and forcing term is estab- 
lished by combining this new representation formula and the method of stationary 
phase. 
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1 Introduction 



This paper is concerned with the initial- and boundary- value problem (IBVP) for the 
Benjamin-Bona-Mahony (BBM) equation in the quarter plane {(x, t),x >0,t> 0}, 

u t + au x + (3uu x - ^u xxt = f, x>0,t>0, 
u(x, 0) = uo(x), x > 0, (1.1) 
u(0,t)=g{t), t>0 

where 7 > 0, a and (5 are real parameters. Our goal is two-fold: first, to establish 
a new representation formula for (1.1) by inverting the operator d t + ad x — ^d xxt \ 
second, to understand if the solution of (1.1) exhibits certain time-periodic behavior 
if the boundary data g and the forcing term / are periodic. 

The study on the large-time periodic behavior is partially motivated by a labo- 
ratory experiment involving water waves generated by a wavemaker mounted at the 
end of a water channel. It is observed that if the wavemaker is oscillated periodically, 
say with a long period T , it appears that in due course, at any fixed station down the 
channel, the wave elevation becomes periodic of period To. Professor Jerry L. Bona 
proposed the problem of establishing this observation as a mathematically exact fact 
about solutions of the suitable model equations for water waves. One goal of this 
paper is to determine if the solution u of (1.1) exhibits eventual periodicity. More 
precisely, we investigate whether the difference 

u(x,t + T ) -u(x,t) (1.2) 

approaches zero as t — > 00 if g is periodic of period T . Since the solution of (1.1) 
grows in time (measured in the norm of Sobolev spaces H k with k > 0)(see [5]), the 
issue of eventual periodicity appears to be extremely difficult. 

The eventual periodicity has previously been studied in several works. In [5], 
Bona and Wu thoroughly investigated the large-time behavior of solutions to the 
BBM equation and the KdV equation including the eventual periodicity. A formula 
representing the solution of the BBM equation was derived through the Laplace trans- 
form with respect to temporal variable t and the eventual periodicity is shown for the 
linearized BBM equation with zero initial data and forcing term. It appears that the 
formula derived there can not be easily extended to include a non-zero initial data 
or a forcing term. Bona, Sun and Zhang in [4] established the eventual periodicity in 
the context of the damped KdV equation 

u t + u x + uu x + u xxx + u = 

with small amplitude boundary data u(0, t) = g(t). They were able to obtain time- 
decaying bounds for solutions of this equation and the eventual periodicity follows as 
a consequence. Through the Laplace transform with respect to the spatial variable 
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x, Shen, Wu and Yuan in a recent work [8] obtained a new solution representation 
formula for the KdV equation and re-established the eventual periodicity of the lin- 
earized KdV equation. In addition, the eventual periodicity of the full KdV equation 
were studied there through extensive numerical experiments. 



In this paper, we first derive a new solution formula for the IBVP 

u t + au x - ~/u xxt = f, x>0,t>0 
u(x, 0) = uo(x), x > 0, 
u(0,t)=g(t), t>0. 

This explicit formula reads 

u(x,t) = g(t)e + / T(x - y,t)u (y) dy 

ft POO 



1.3) 



o 



+ 

'o Jo 

where T and $ are given by 



$(x-y,t-T) 



a 



f(y,r) + —g(r)e ^ 
V7 



dydr, (1.4) 



F(x,t) 



f 

J — c 



e !+^« 2 a4, 



/OO 1 ^ 



T should be understood as a distribution. To obtain (1.4), we consider both the even 
and odd extensions of (1.3) to the whole spatial line. By taking the Fourier transforms 
of these extensions and solving the resulting equations simultaneously, we are able to 
represent 



sin(a;£) u(x, t) dx and 



cos(x^) u(x, t) dx 



1.5) 



in terms of /, Uq and g. (1.4) is then established by taking the inverse Fourier 
transform of the quantities in (1.5). Corollaries of (1.4) include explicit solution 
formulas of (1.3) with a = and 7 > and of (1.3) with 7 = 0. More details can be 
found in the second section. 

In [1] and [3], the IBVP (1.1) has been recast as an integral equation through the 
inversion of the operator d t — "yd xxt , 



u(x,t) =u (x) + g(t)e ^ + 



t POO 



JO 



K(x,y) (au+ -f3u 2 ) {y,r)dydr, (1.6) 



where the kernel function K (x, y) is given by 

1 r 



K(x,y) 



2 7 



x + y 



e V7 4- sgn(x — y)e 



:i-7) 
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While this representation is handy in dealing with the well-posedness issue, we find 
it inconvenient in studying the eventual periodicity of (1.1) due to the inclusion of 
the linear term on the right. The new representation (1.4) allows us to show that any 
solution of (1.3) is eventual periodic if / and g are periodic of the same period. The 
precise statement and its proof are presented in the third section. 

We mention that there is adequate theory of well-posedness and regularity on the 
IB VP (1.1). The following theorem of Bona and Luo [2] serves our purpose well. In 
the following theorem, we write R + = [0, oo) and C% is exactly like C k except that 
the functions and its first k derivatives are required to be bounded. 

Theorem 1.1 Let I = [0, T] if T is positive or I = [0, oo) if T — oo. Assume that 
g G C 1 (J) and uq G Cf(R + ) fli/ 2 (R + ). Then (1.1) is globally well-posed in the sense 
that there is a unique classical solution u G C l (I, C£°(R + )) n C(I; H 2 (R + )) which 
depends continuously on g G C 1 (J) and uq G C|(R + ) fl H 2 (R + ). 

A more recent work [1] reduces the regularity assumptions to g G C(J) and uq G 
C^(R + ) while (1.1) still has a unique global solution in a slightly weak sense. We 
shall not attempt to optimize these regularity assumptions in this paper. The rest of 
this paper is divided into two sections and two appendices. 



2 The inversion of the operator d t + ad x — jd xx t 

This section explicitly solves the IBVP of the linearized BBM equation 

u t + au x - ~fu xxt = f, x > 0, t > 0, 

< u(x, 0) = u (x), x > 0, (2.1) 

k u(0,t) = g(t), f >0, 

where 7 > and a are real parameters. This amounts to inverting the operator 
d t + ad x — jd xxt for the quarter-plane problem. Without loss of generality, we assume 
9(0) =0. 

In the case when 7 > 0, we consider a new dependent variable 

w(x,t) = u(x,t) -g(t)e~^, (2.2) 
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which satisfies the equations 



w t + aw x - 7w xxt = f, x > 0, t > 
w(x, 0) = u (x), x > 0, 
k w(0,t) = 0, t > 0. 



(2.3) 



where 



/(x,t) =/( a: ,i) + - F ^)e ^- 



(2.4) 



The IBVP (2.3) has its boundary data equal to zero and is solved through odd and 
even extensions to the whole spatial line. The solution u of (2.1) is then obtained by 
(2.2). 

Theorem 2.1 Let I = [0, T] if T is positive or I = [0, oo) if T = oo. Let u e 
Cf(R+) n # 2 (R+) ; G C^J) and f e C(/,L 2 (R+)). W^/w^ loss of generality, 
assume tto(0) = g(0) = 0. T/ien the unique classical solution of (2.1) can be written 
as 



u(x,t) = g(t)e + / T(x -y,t)u {y)dy 

Jo 

ft poo 

I / <!'(./• //./ v) 



'o Jo 

where T and $ are given by 



a 



f(y,r) + —g(r)e ^ 

V7 



r(x,t) = 



ix^—i 



i+W df, 



■oo 



ix£—i 



T m ^.^) should be understood as a distribution. 



dj/dr, (2.5) 

(2.6) 
(2.7) 



To gain an initial understanding of the formula in this theorem, we consider two 
special cases. The first is a = and / = 0. When a — 0, 



r(x,f) = / e ia *dZ = 6(x), 

J — oo 

where 5 denotes the Dirac delta. Therefore, for x > 0, 

POO 

/ T(x-y,t)u (y)dy = u (x). 
Jo 
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Corollary 2.2 The solution of (2.1) with a = and f = is given by 

x 

u(x,t) — u (x) + g(t)e Vt\ 

The second special case is when 7 = and / = 0. Although 7 = in not allowed 
in Theorem 2.1, the solution formula for this case can still be obtained similarly as 
(2.5). Instead of (2.2), one considers 



which solves 



w(x, t) = u(x, t) — g(t)e x 

w t + aw x = (ag(t) - g'(t))e- x , x > 0, t > 0, 

w(x, 0) = u (x), x > 0, 

k w(0,t) = 0, t > 0. 
Then, as in Theorem 2.1, the solution for this special case is 

PCX) 

u(x,t) = g(t)e' x + / T(x-y,t)u (y)dy 

Jo 

+ 11 $(x-y,t-T)[ag(r)e~ y - g\r)e- y ] dydr. (2.8) 
Jo Jo 

This representation allows us to extract a simple solution formula for the IBVP (2.1) 
with 7 = and / = 0. 



Corollary 2.3 The solution of the IBVP (2.1) with 7 = and f = is given by 

u(x, t) = 



u (x — at), if x > at, 
g (t j , if x < at. 

Proof of Corollary 2.3. When 7 = 0, 



/oo 
■00 

If x — at > 0, then 

POO POO 

/ T(x -y,t)u (y)dy = 5(x - at - y)u (y) dy = u (x - at) (2.9) 
Jo Jo 
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and 

"t poo 



pi rco 

/ / $(x -y,t -t) [ag(r)e~ y - g\r)e" yl \ dydr 
Jo Jo 

pt POO 

= / [a^(r) — g'(r)] / 5(x — at + ar — y) e~ y dy dr 
Jo Jo 

= f \ag{r) - g\r)\ e -(*-<**") rfr 
./o 

/ «^(r)e- aT rfr- / g'^e'^dr 
Jo Jo 



—x+oct 



= -9(t)e- x . (2.10) 
Inserting (2.9) and (2.10) in (2.8) yields 

u(x, t) = u (x — at) if x — at > 0. 

If x — at < 0, then 

/•oo 

/ 5(x - at - y)u (y) dy = 0. 
Jo 



and 

0, otherwise. 



y 8(x — at + ar — y) e y rfy = | 



Thus, 

pt /»oo 

/ / $(x — y, t — r) [ag(r)e~ v — g'(r)e~ y ] dydr 
Jo Jo 

= f e -^-^r) N(r) _ g>{r)] dT 

J t—a~ 1 x 

= -g(t)e- x + g{t-a- 1 x). 

Therefore, for x < at, 

u(x, t) = g(t — a~ x x). 
This completes the proof of Corollary 2.3. 

Proof of Theorem 2.1. The major idea is to extend the equation of (2.3) from the 
half line {x : x > 0} to the whole line xGRso that the method of Fourier transforms 
can be employed. Both odd and even extensions will be considered. The rest of the 
proof is divided into four major steps. 



Step 1. Odd extension. Denote by W the odd extension of w in x, namely 
W(x,t) = 



w(x, t) if x > 0, 

— w(— x, t) if x < 0. 



If Wo and F are the odd extensions of wo and /, respectively, then W solves the 
following initial-value problem 



(2.11) 



W t + a sgn(x)W x - -fW xxt = F, xeR,t>0 
W(x,0)=W o (x), xeR. 
Let W denote the Fourier transform of W , namely 

/oo 
e-* x W{x,t)dx. 
-oo 

After applying a basic property of the Fourier transform, we obtain 

(l + 7f)9 t %i) + «^Wa;)^) = %,t). (2.12) 
According to the definition of W, we have 

roc poo 

W(£,t)= (e~ ix t - e ix t)w(x,t)dx = -2i sm(x£)w(x,t)dx. (2.13) 
Jo Jo 

In addition, 

/0 poo 
e-^ x sgn(x)W x dx+ / e-^ x sgn(x)W x dx 
oo Jo 

/0 poo 
e~^ x (-w x (-x, t))dx + / e~^ x w x (x,t)dx. 
oo Jo 

Making the substitution y = —x and integrating by parts yield 

poo POO 

F(sgn(x)W x ){{;,t) = - e^w y (y,t)dy+ e^ x w x (x,t)dx. 

Jo Jo 

POO 

= iU (e* 3 * + e~ i3 *)w(x,t)dx 
Jo 

POO 

= 2i£ I cos(x£) w(x,t)dx. (2.14) 
Jo 

We can also write F in terms of / as 

^00 

F(£,t) = -2i / sm(xOf(x,t)dx (2.15) 
Jo 

Inserting (2. 13), (2.14) and (2.15) in (2.12), we obtain 

d t X(U)-mY(U) = (2-16) 
9 



where 

poo poo 

X(£,t) = / sin(x^) w(x, t) dx, Y(£,t) = / cos(x^) w(x, t) dx, 
Jo Jo 

1 f°° 

P(Q = 1 + ^ 2 and fei(£,t) = - - — 2 J sm(x£) f(x,t)dx. 

Step 2. Even extension. Denote by V(x,t) the even extension of w, namely 

, / w(x,t) if x > 0, 

«;(-*,*) ifx<0. 

It can be easily verified that V satisfies 

V t + a sgn(x)V x - ^V xxt = H, x G R, t > 

(2.17) 

V(x,0) = V o (x), xeR, 

where H and Vo are the even extensions of / and wo, respectively. As in Step 1, we 
have 

^ r°° 

V(£,t) = 2 cos(x£) w(x,t)dx, (2.18) 
Jo 

and 

poo 

H(£,t) = 2 cos{xOf(x,t)dx (2.19) 
Jo 

Furthermore, 

e-^sgn(x)\4rfx+ / e-* x sgn(x)V x dx 

-oo 

/O /"DC 
e~^ x w x (-x,t)dx + / e~^ x w x (x,t)dx. 
-oo Vo 

Making the substitution y = — x in the first integral and integrating by parts leads to 

/■oo /*oo 

.F(sgn(x)\4)(£,t) = / e* v w y (y,t)dy+ e-* x w x (x,t) dx 

Jo Jo 

poo poo 

= ~i£ / e l ^ x w(x, t) dx + / e~^ x w(x,t) dx 
Jo Jo 

poo 

= 2f / sin(x£) w{x, t) dx. (2.20) 
Vo 

Taking the Fourier transform of (2.17) and applying (2. 18), (2.19) and (2.20), we 
obtain 

d t Y(Z,t)+(3X(Z,t) = h 2 (Z,t), (2.21) 
10 



where 



1 f°° 

h 2 (£,t) = i + y cos(x£) /(x, da;. 



Step 5. SoZuira^ for X(£,t) and Y(£,t). Solving the system of (2.16) and (2.21), 
we find 

roc pt poo 

X(£,t)= sin(x£ + (3t) u (x)dx + - — / / sin(^ + (3{t - r)) f(x, r)dxdr, 

Jo 1 + 7s Jo Jo 

/• \ j pt POO ^ 

Y(£, t)= cos(x£ + pt) u (x)dx H / / cos(^ + (3(t - r)) /(x, r)dxdr. 

Jo 1 - 7<r ./o ./o 

where / is defined in (2.4). We leave the details of derivation in Appendix A. 

Step 4- Finding w(x, t) through the inverse Fourier transform. To find w(x, t), we 
first notice that 

POO 

/ e- i *w(y,t)dy = Y(Z,t)-iX(U)- 
Jo 

Applying the formulas in the previous step, we have 

POO POO 

/ e-^w(y,t)dy = e' if3t / e-^u (y)dy 
Jo Jo 

i pt POO 

+TTW 2 Jo e ~ m ' T) Jo e ^tt y ^ dydT - ( 2 - 22 ) 

We shall now establish a theorem that allows us to obtain w(x, t) by taking the inverse 
Fourier transform of (2.22). 

Theorem 2.4 Fix t > 0. Ifu(x,t) G L 2 (R + ) ; then, for any x > 0, 

/OO POO 
jxti i e - iy t U (y,t)dyd£ = u(x,t). (2.23) 
oo Jo 

Proof of Theorem 2.4- Recall that if g e (£) = exp (— e7r£ 2 ), then 

^ 1 / TTX 2 

g e (x) = ^exp ^ — 

In addition, for any / G L P (R) with 1 < p < oo, 

&*/ - / in LP(R). 
These basic facts can be found in Lieb and Loss [6]. Now, consider 

»oo poo 



/OO /' 
■oo JO 



e^(0 / e-^(y,0rfy^- 
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According to Lemma 2.5 below, if u(y,t) G L 2 (R + ), then 

POO 

P(u)(U)= / e-^u(y,t)dyeL 2 (R). 



Since g e (Q — ► 1 as e — > 0, the dominated convergence theorem implies 

&P(u) -> P(u) in L 2 (R). 

Then, 

/OO /»00 /»00 /»00 

e ia *&(0 / e-^(y,0rfy^ = / <y,t) / e ~^%(0 ^2/ 
-oo JO JO J-oo 

POO 

= / u(y,t)g e (y-x)dy. (2.24) 
./o 

As in the proof of Theorem 2.16 of [6], we can prove 

roc 

/ u(y,t)g e (y - x)dy -> in L 2 (R+). 

Jo 

Letting e -> in (2.24) yields (2.23). This proves Theorem 2.4. 
Lemma 2.5 Ifu(y,t)e L 2 (R + ), then 

POO 

P(u)(U)= / e-^u(y,t)dyeL 2 (R). 
Jo 

Proof of Lemma 2.5. For any e > 0, 

/oo 
P(u)(U)P(u)(UMOdti 
-oo 

oo poo poo 

g e (0 / e-^u{x,t)dx / j*u(y,t)dyd£ 
Jo Jo 



oo 

noo POO 
u(x,t)u(y,t) / e-* x -Mg e (Z)dtdxdy 
J-oo 

POO POO 

= / u(x,t) / g e (x - y)u(y,t)dydx 
Jo Jo 

As e -> 0, 

/•oo 

/ de(x -y)u(y,t)dy -> w(x,t) in L 2 (R + ) 
./o 

and u(x,t) £ L 2 (R + ) implies J*^ |P(w)| 2 (£, t)g e (£) d^ is bounded uniformly. Since 

= exp(-e7r|£| 2 ) -> 1 as e -> 0, 

12 



we obtain by applying the monotone convergence theorem that 



\P(u)\ 2 (£,t)dS 



\u(x, t)\ 2 dx. 



This proves Lemma 2.5. 



Taking the inverse Fourier transform (denoted by T v ) of (2.22) and applying 
Theorem 2.4 and the basic property T~ x {fg) = ^ r_1 (/) * .F _1 (<7), we obtain 



w(x, t) 



Y(x - y,t)u (y) dy + 



t POO 



JO 



$(x -y,t- r)f(y, r) dydr, 



where, noticing f3 = jf^, 



ix£—i . \» t 



r(x,t) = / e^^'dZ, 



■oo 

OO i 



1 + 7^ 



■ A ■ Ctlr . 

-e 1+ t« 2 a£. 



Therefore, by (2.3), 

— x 

u(x,t) = w(x,t)+g(t)e ^ 

poo 

= g(t)e~si + F(x-y,t)u (y)dy 
Jo 



pt poo 

+ / / $(x-y,t-r) 
Jo Jo 

This completes the proof of Theorem 2.1. 



cx y 
f(y,r) + —g(r) e'^ 
V7 



dydr. 



3 Eventual periodicity 

This section studies the eventual periodicity of solutions to the IB VP for the linearized 
BBM equation 

u t + au x - ju xxt = f(x, t), x > 0, t > 0, 
u(x,0) = u (x), x>0, (3.1) 
u(0,t) = g(t), t>0 

where g and / are assumed to be periodic in t. We prove the following theorem. 
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Theorem 3.1 Let 7 > and a be real parameters. Let I = [0, T] if T is positive or 
I = [0,oo) ifT = 00. Letu E C 6 2 (R+) andu' E L 1 (R+). Let g E C{I) with g(0) = 
and f E C(I, L 1 (R + )). Assume g and f are periodic of period T in t, namely, for 
all t > 0, 

g(t + T )=g(t) and f(x,t + T ) = f(t). (3.2) 
Then, for any fixed x > 0, the solution u of (3.1) satisfies 



lim (u(x, t + T ) - u(x, t)) = 0. 



t— >oo 



(3.3) 



That is, u is eventually periodic of period T . 



Remark 3.2 When 7 = 0, the eventual periodicity is easily obtained from the explicit 
formula in Corollary 2.3. 



Proof of Theorem 3.1. Consider the new function 

v(x, t) = u(x, t) — Uq(x), 

which satisfies 

v t + av x - jv xxt = f(x,t) -au' (x), x>0,t>0, 
v(x,0) = 0, x>0, 
v(0,t) = g(t), t>0 

Applying the explicit formula to (3.4) gives 



(3.4) 



v(x,t) = g(t)e 

pt poo 

+ / / $(x-y,t-T) 
Jo Jo 



a 



f(y,r)-au' (y) + —g{T)e ^ 



dydr, 



where $ is given by (2.7). Noticing the conditions in (3.2), we obtain after making a 
substitution, 



u(x, t + T ) — u(x, t) = v(x, t + T ) — v(x, t) 

»t+T poo 



pt+i poo 

/ / $>(x-y,t + T -T) 
Jo Jo 



a 



f(y,r)-au' (y) + —g(r)e ^ 



dydr 



t POO 



<$>(x-y,t-r) 



JO 

poo 



a 



f(y,r) - au' Q (y) + —g{r)e v=r 



T JO 



<f>(x-y,t-r) 



a 



f(y,T)-au' (y) + —g(T)e ^ 



dydr 

dydr. (3.5) 



To show (3.3), we first show that, for any e > 0, there is K > such that 

COO J 



Mx,t)\ 



ix£—i - t 7/ . 

e !+^« 2 d£ 



< e when t > K. 
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First, we choose M = M(e) > such that 



/— M j poo | 



Next, we apply the method of stationary phase to show the following asymptotics. 
Proposition 3.3 For any fixed M > 0, 



/M 
■M 



l + T^ 7 di = O 



1 

71 



- M i + ie 

as t — > oo. 77ms large-time asymptotics is uniform in x e R + . 

We remark that the result of this proposition does not necessarily hold for 



7 — c 



ix^—i 



at 



1+7*2 



because one of the conditions in the method of stationary phase is violated. More 
details on this point will be provided in Appendix B. The following estimate is a 
special consequence of this proposition. 

Corollary 3.4 There exists K = K(M) such that 



M 



1 



< 



- M i + ie 

whenever t > K . 

We resume the proof of Theorem 3.1. It then follows from (3.5) that, for t > K, 



\u(x, t + T) — u(x, t)\ 

r roo 



< e 



poo pu 

\f(y,r)\dydr + aT / \u' (y)\dy + a / g(r)dr 

JO J-T 



f -T JO 

This completes the proof of Theorem 3.1. 

Proof of Proposition 3.3. We apply the method of stationary phase. To do so, we 
split the integral into four parts 



/M 
-M 



1 



i+^'df = h+I 2 + h + h 



+ 



M 



+ 



+ 



M 
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Note that ±1/^/7 are the zero points of the derivative of = a£/(l+7£ 2 ) and p'(0 
is nonzero for £ in each of these intervals. It suffices to consider I 2 and I4. Without 
loss of generality, we assume a > 0. Direct applications of the method of stationary 
phase concludes that 

^ 6 4 6V " (3.6) 

(3.7) 



For readers' convenience, this method is recalled in Appendix B and the details leading 
to these estimates can also be found there. 
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Appendix A 

This appendix provides the details of solving the system of ODEs (2.16) and (2.21). 
Consider the general nonhomogeneous linear systems 

^x(t) = P(t)x(t) + g(t), x(t ) = x 

where x e R™, P £ R nxn and g e R n . By variation of parameters, its solution can 
be written as t 

x(f) = ^(t)^- 1 (t )x + ^(t) / ^- 1 (s)g(s)ds (A.l) 

J to 

where ^ denotes a fundamental matrix of the homogeneous system 

|x(t) = P(f)x(f). 
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Since the system of equations we are solving can be written as 



p 
-p 



x(U) 



+ 



and a fundamental matrix of the corresponding homogenous system is given by 



*(*) 

we apply (A.l) to obtain that 
= tt(t)tf _1 (0) 



i/3t 



X&O) 



Jo 



dr. (A.2) 



Inserting 
and 



e ~i/3t _i e -i@t 







1 





sin(x^)u (x) dx 

i 

cos(x£)uq(x) dx 



in (A.2), we find after some simplification 



1 



sm(x£)f(x, t) dx 



i f°° 

2 / cos(x£) f(x,t)dx 



t POO 



f°° 1 

X(£,t)= sm(x£ + pt)u (x)dx + -—- 

Jo 1 + 7? Jo Jo 



sin(a;£ + P(t — r)) f(x, r)dxdr, 



Y(U) 



f 

Jo 



cos(a;£ + pt) u (x)dx 



1 



l+7^ 2 



t POO 



cos(x£ + P(t — r)) f(x, r)dxdr. 



o jo 



Appendix B 

This appendix offers an expanded commentary on the asymptotic analysis of the 
oscillatory integrals arising in section 3. These analysis relies upon standard results 
in the theory of stationary phase, e.q. Theorem 13.1 in F. Olver's book [7]. For 
readers' convenience, we first recall this theory here. 



Suppose that in the integral 



I(t) 



q(y) dy 
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the limits a and b are independent of t, a being finite and b(> a) finite or infinite. 
The functions p(y) and q(y) are independent of t, p(y) being real and q(y) either real 
or complex. We also assume that the only point at which p'(y) vanishes is a. Without 
loss of generality, both t and p\y) are taken to be positive; cases in which one of them 
is negative can be handled by changing the sign of i throughout. We require 

(i) In (a, b), the functions p'(y) and q(y) are continuous, p'(y) > 0, and p"(y) and 
q'(y) have at most a finite number of discontinuities and infinities. 

(ii) As y — > a+, 

p(y) - p(a) ~ P(y - a)", q(y) ~ Q(y - a)*" 1 , (B.l) 

the first of these relations being differentiable. Here P, /i and A are positive 
constants, and Q is a real or complex constant. 



(iii) For each e G (0, b — a), 



/g(y)V 



< oo. 



(iv) As t — > 6—, the limit of q(y)/p'(y) is finite, and this limit is zero if = oo. 

With these conditions, the nature of asymptotic approximation to I(t) for large t 
depends on the sign of A — fj,. In the case A < /i, we have the following theorem. 

Theorem B.l In addition to the above conditions, assume that A < /i, the first of 
(B.l) is twice differentiable, and the second of (B.l) is differentiable, then 

Q„/A\ e itp{a) 



i(t) ~ e AW(2M)^ r (:M 



fj, KfiJ (Pt) A /M 

as t — > oo. 



We now provide the details leading to (3.6). It suffices to check the conditions of 
Theorem B.l. Setting 

:, 6 = 0, p(0 = _*- and q(0 



we have 

(i) p, p', p", q and q' are all continuous in (—1/^/7,0), and p'(0 > 0- 



18 



(n) as e - -3=+, 



p(0-p(-^ 



2 ^ + ^f) ' 



5 



1 

-e . 
2 



That is, P = ^p-, n = 2,Q = \e and A = 1. 
(iii) For any fixed e > 0, V_j_ +t (q/p') < oo. In fact, 



a(l- 7 £ 2 ) ' ^ £ 'Vi y_^ +e 



d£ < oo. 



(iv) As £ -> 0-, g/p' -> 1/a. 
Theorem B.l then implies 
,o j 



e s i+7« 2 rf£ ~ e * _ e vi r - 



i+ 7 e 



2; (a/2^) 1/2 



The estimate (3.7) also follows from Theorem B.l. The conditions can be similarly 
checked for this integral. In fact, for 



a = — , b=M, p(0 



a/7 



l + 7 £ 5 



and g(0 



1+7^ 



we have 



1 i^x 

~ — e 
2 



It is also easy to verify that Vj_ +e m(q/p') < 00 f° r an y fixed e > 0. In fact, 

s/7 ' 



P 



a(l — 7 £ 2 ) ' 



p 



M 



q\> 
p' 



d!; < 00. 



(B.2) 



In addition, as £ — > M, q/p' tends to a finite limit. It then follows from Theorem B.l 
that 



e a£ ~ 



When M = 00, (B.2) can not be verified and Theorem B.l does not apply to the 
integral on (—00,00). 
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